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Abstract 

We calculate the supercurrent through a Josephson junction consisting of a phase-coherent metal 
particle (quantum dot), weakly coupled to two superconductors. The classical motion in the quan- 
tum dot is assumed to be chaotic on time scales greater than the ergodic time r erg , which itself is 
much smaller than the mean dwell time Td we ii- The excitation spectrum of the Josephson junction 
has a gap -E ga p; which can be less than the gap A in the bulk superconductors. The average 
supercurrent is computed in the ergodic regime r er g <C h/A, using random-matrix theory, and 
in the non-ergodic regime r er g 3> h/A, using a semiclassical relation between the supercurrent 
and dwell-time distribution. In contrast to conventional Josephson junctions, raising the tem- 
perature above the excitation gap does not necessarily lead to an exponential suppression of the 
supercurrent. Instead, we find a temperature regime between E gSuP and A where the supercurrent 
decreases logarithmically with temperature. This anomalously weak temperature dependence is 
caused by long-range correlations in the excitation spectrum, which extend over an energy range 
h/T evg greater than E g£Lp ~ ?i/ T dweii- A similar logarithmic temperature dependence of the super- 
current was discovered by Aslamazov, Larkin, and Ovchinnikov, in a Josephson junction consisting 
of a disordered metal between two tunnel barriers. 
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I. INTRODUCTION 



The dissipationless flow of a current through a superconductor-normal-metal- 
superconductor (SNS) junction is a fundamental demonstration of the "proximity effect": a 
normal metal borrows characteristic properties from a nearby superconductor. The energy 
gap A in the bulk induces a suppression of the density of states inside the normal metal 
near the Fermi level, depending on the phase difference between the superconductors. The 
resulting ^-dependence of the free energy F implies the flow of a current / = (2e/h)dF/d<j) in 
equilibrium. In contrast to the original Josephson effect in tunnel junctions, the separation 
of the superconductors in an SNS junction can be much greater than the superconducting 
coherence length. Recent experiments on mesoscopic Josephson junctions [1, 2, 3, 4, 5, 6] 
have revived theoretical interest in this subject [7, 8, 9], which goes back to work by Kulik 
[10] and Aslamasov, Larkin, and Ovchinnikov [11]. (For more references, see the review 12.) 

In this paper, we consider the case that the normal region consists of a chaotic quantum 
dot. A quantum dot is a small metal particle, within which the motion is phase coherent, 
weakly coupled to the superconductors by means of point contacts. We assume that the 
classical dynamics in the quantum dot is chaotic on time scales longer than the time r crg 
needed for ergodic exploration of the phase space of the quantum dot. (In order of mag- 
nitude, r crg ~ L/vp for a quantum dot of size L without impurities, v-f being the Fermi 
velocity.) On energy scales smaller than ^/r erg , the spectral statistics of a chaotic quantum 
dot is described by random-matrix theory [13, 14]. On larger energy scales, the non-ergodic 
dynamics on time scales below r erg becomes dominant [15]. The condition of weak coupling 
means that the mean dwell time Td we ii in the quantum dot is much greater than r crg . (The 
ratio Tdwcii/Tcrg is of the order of the ratio of the total surface area to the area of the point 
contacts.) Although Josephson junctions are commonly known as "weak links" [16], we will 
refer to junctions where Tdweii — fcrg as "strongly coupled" junctions, to distinguish them 
from the weakly coupled junctions (rdwcii 3> T crg ) considered here. 

A weakly coupled SNS junction consisting of a dirty normal metal separated from the two 
superconductors by high tunnel barriers was studied in the original paper by Aslamasov, 
Larkin, and Ovchinnikov [11]. Their theory was restricted to the high-temperature regime 
kT ^> ^/r dwell . In contrast to strongly coupled Josephson junctions, where the supercurrent 
is suppressed exponentially for kT > h/r dwc n [9, 17], it was found that the supercurrent 
depends logarithmically on temperature for kT < h/r erg , while exponential suppression 
only sets in when kT > h/r crg . In the present paper we find a qualitatively similar high- 
temperature behavior of the supercurrent in the case that the weak coupling is ensured by 
point contacts rather than tunnel barriers. In addition, we are able to go down to zero 
temperature, where we find that the supercurrent acquires a logarithmic dependence on 
the minimum of Td W eii/T erg and Td we iiA/fr, over and above the conventional dependence on 
min(^/rd W cii, A) known from strongly-coupled Josephson junctions. 

Our paper builds on earlier work with Melsen and Frahm [18], where we computed the 
density of states of a chaotic quantum dot which is weakly coupled to a superconductor, 
and found a gap at the Fermi level of width E gSup ~ h/rd we \\ (provided Ti/r^w <C A). 
Although the supercurrent can be expressed as an integral over the density of states, direct 
application of the results of Ref. 18 is not possible, as they were only derived under the 
condition e % ^/rd we ii. We have it found it necessary to relax this restriction, because the 
sensitivity of the excitation spectrum to the phase difference of the superconductors (which 
determines the supercurrent) extends up to Ti/ r crg . It is because of this long-range sensitivity, 
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that the temperature scale for the exponential suppression of the supercurrent is set by r crg 
and not by r dwell . 

This work is organized as follows: In Sec. II we review the scattering theory of the 
Josephson effect [19, 20], on which our calculation is based. We distinguish two regimes, 
the ergodic and non-ergodic regime, depending on the relative magnitude of A and h/r erg . 
In Sees. Ill and IV, we consider the ergodic regime r crg <^h/A, where we can use random- 
matrix theory. The non-ergodic regime r erg > h/A is treated in Sec. V, using a semiclassical 
relation between the supercurrent and the dwell-time distribution. We conclude with an 
overview of our results in Sec. VI. 



II. SCATTERING MATRIX FORMULA FOR THE SUPERCURRENT 

We consider the SNS junction sketched in Fig. 1. The two superconductors Si and S2 
have order parameters Ae^ 1 , Ae*^ 2 , with phase difference = 0i — 02. The contacts to 
the normal metal N have Ni, N 2 propagating modes at the Fermi energy E-p. We denote 
N = Ni + N 2 . Elastic scattering by the normal metal at energy E = E F + e is characterized 
by an N x N unitary matrix S(e). Excitations in N with energy e > consist of electrons 
(occupied states lying e above the Fermi level) and holes (empty states lying e below the 
Fermi level). Their scattering matrix S-^(e) has dimension 2N x 2N, with the block structure 

~,,(S(e) \_(S"S*\ m 

6 N (£j- ^ 5 * ( _ £) j = y S he S hh y (1) 

The off-diagonal blocks are zero, because the normal metal does not scatter electrons into 
holes. 

The supercurrent couples electron and hole excitations through the mechanism of Andreev 
reflection [21]: An electron approaching the NS interface from the normal side at e < A is 
reflected as a hole, and vice versa. The scattering matrix Sa(z) for Andreev reflection is 
given by (assuming A <C E F ) 

( e i$ \ 

S A (e) = a(e) , $ , (2a) 





a(e) = e" iarccos(£/A) = e/A - iy/l - s 2 /A 2 . (2b) 

Here $ is a diagonal matrix with elements $jj = 4>i = 0/2 for 1 < j < Ni and = 2 = 
—0/2 for JVi + 1 < j < Ax + N 2 . The matrix Sa has been defined for e < A. Its definition 
may be extended to e > A, when ct(e) = e/A — \/e 2 — A 2 . Notice that Sa is no longer 
unitary for e > A. 

The matrices <Sn and Sa determine the excitation spectrum of the Josephson junction for 
all e, in the following way [19]. For e < A, the spectrum is discrete. The density of states 
p(e) = J2n^( e ~ e n) consists of delta functions at the solutions of the equation 

det [1 - S A (e n )S N (e n )] = 0. (3) 

For e > A, the spectrum is continuous, with density of states 



p(e) = --Im-^- 
7r de 



lndet [1 — 5a(£)5n(£)] — ^ lndet Sa (e)S-n(e) 



(4) 
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FIG. 1: Schematic drawing of a chaotic Josephson junction. 



If e is replaced by e + i0 + , Eq. (4) describes both the continuous and the discrete spectrum 
]. The excitation spectrum determines the free energy F of the Josephson junction, 



poo 

F = -2kT / dep(e) In [2 cosh(e/2kT)] + [^-independent terms], 
Jo 



(5) 



where T is the temperature. The ^-independent terms include [23] a spatial integral over 
|A(r*)| 2 , which does not depend on in the step-function model for the pair potential A(r). 
Only the ^-dependent terms contribute to the supercurrent 



_ 2edF_ 

h d(j> 



(6) 



We use the analyticity of <Sa and <Sn in the upper half of the complex e-plane to rewrite 
the expression for the supercurrent in a more convenient form. Under a change e — > —e, the 
determinant det[l — Sa(z + ?0 + )5n(£ + ^0 + )] goes over into its complex conjugate. Combi- 
nation of Eqs. (4)-(6), and extension of the e-integration from — oo to oo, results in 



6 -2krl 



oo+i0+ 



d 



, .. <fe In [2 cosh(e/2fcT)] — lndet[l - 5 A (e)5 N (e)]. 

irch d<p i_ 00+i o+ de 



(7) 



We now perform a partial integration and close the integration contour in the upper half 
of the complex plane. The integrand has poles at the Matsubara frequencies iuj n = (2n + 
l)mkT. Summing over the residues one finds 



2e d °° 
/ = -42kT— Vlndetfl - S A {tu n )S N {tuj n )}. 
n d(p ' 



(8) 



n=0 



Eq. (8) is the starting point for our evaluation of the average supercurrent through a chaotic 
Josephson junction. 



III. SUPERCURRENT THROUGH A CHAOTIC JOSEPHSON JUNCTION 

We consider the case that the normal region has a chaotic classical dynamics on time 
scales greater than the ergodic time r crg . In this section, we assume that r crg h/A, so 
that we may use random-matrix theory to evaluate the ensemble average of the supercurrent. 
We postpone to Sec. V a discussion of the regime r erg ^ h/A, in which the non-ergodic 
dynamics on time scales shorter than r erg starts to play a role. We assume that the normal 
metal is weakly coupled to the superconductors, so that the mean dwell time Td we ii ^> r erg . 
No assumption is made regarding the relative magnitudes of r dwc n and h/A. 
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We use a relationship between the scattering matrix S of the normal metal and its 
Hamiltonian H [24, 25], 

S(e) = 1 - 2mW\e - H + iirWW^W. (9) 

The Hamiltonian H (representing the isolated normal metal region) is taken from the Gaus- 
sian ensemble of random-matrix theory [26] , 

P(H) oc exp ^-^MA _2 tr H 2 ^j , (10) 

where M is the dimension of H (taken to infinity at the end) and A is a parameter that 
determines the average level spacing 5 = \tt/2M of the excitation spectrum in the normal 
region. (This spacing 5 is half the level spacing of H, because it combines electron and hole 
levels together.) The matrix H is real and symmetric. The coupling matrix W is an M x N 
matrix with elements[27, 28] 



= U mn {2M5) 1 ' 2 (2T~ l - 1 - 2r-Vi-r„) V2 . (11) 

Here Y n is the transmission probability of mode n in the contacts to the superconductor. 
For ballistic contacts, T n = 1, while T n <ti 1 for tunneling contacts. 

We now substitute Eq. (9) for S into Eq. (1) for and then substitute 5n into Eq. (8) 
for the supercurrent. Using also Eq. (2) for Sa, we find after some straightforward matrix 
algebra that 



oo 



/ = -^2fcT— ^lndet[iw n -W + W(iw n )], (12) 

™ n=0 

where we have introduced the 2M x 2M matrices 

h =( H o-h)> ^ 

Me) ~ nA ( ^/ A ) WWT We '* WT \ fl3b) 

The matrix 7i — W(e) is the effective Hamiltonian of Refs. 18, 29 (where the regime c<C A 
was considered, in which the e-dependence of W(e) can be neglected). 
We define the 2M x 2M Green function 

0(e) = [e-n + W(e)}-\ (14) 

which determines the density of states according to 

p(e) = -ir-HmtrQie + iO) (l + . (15) 

Eq. (15) is equivalent to Eq. (4). The expression for the supercurrent in terms of Q(e) is 

2 d °° 
I = -^2kT—y^\ndetg(iu n ) 

n dcp 

2e °° d 
= -^-2kTj2^G(ico n )j-W(iu; n ). (16) 

n=0 ^ 
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The average supercurrent follows from the average Green function (Q(e)), since W is a 
fixed matrix. The average over the random Hamiltonian H (determining Q) is done with the 
help of the diagrammatic technique of Refs. 30, 31. We consider the regime M, N, \e\/5 ^> 1 
in which only planar diagrams need to be considered. Resummation of these diagrams leads 
to a self-consistency equation which is similar to Pastur's equation [32], 

(0(e)) =[e + W(e) - (X 2 /M)V(e) ® 1 M ] 1 . (17) 

The symbol <8> indicates the direct product between the M x M unit matrix 1 M and the 
2x2 matrix 

/ (trS-) _( tr 0<*)\ 



-(trg he ) (trg hh ) 
We seek the solution of Eqs. (17)— (18) which satisfies 

eG(e) -> 1 2M if \e\ > A. (19) 

It is convenient to define a self energy 

trQ ee trG eh 



f s ee 






1 S fee 




: M \ 



M \ ^G he tiQ hh 1 ^ 

Eqs. (17)-(20) contain a closed set of equations from which (S) can be determined. We are 
interested in the limit M — > oo, A — > oo, keeping N and e/5 = 2eM/\n fixed. In this limit, 
the equations for (£) become 

(S ee ) = (S M ), (S e/l )(S /ie ) - (S ee ) 2 = 1, (21a) 

TV 

^(•E eh ) +J2 K i {z(^ eh ) + Ae^(S ee )) = 0, (21b) 

N 

^£( S ^) + Y^Kj (e{Y, he ) + Ae-^(S ee )) = 0. (21c) 

3=1 

The function Kj(e) is defined through 

T j /K J = (4 - 2r j )v / A 2 - e 2 + Fj (Ae _i *" (S e/l ) + Ae 1 ^' (S he ) + 2e(S ee )) . (21d) 

(We substituted Eq. (11) for the matrix W.) The boundary condition (19) becomes ineffec- 
tive in the limit A — > oo. Instead, we seek the solution of Eq. (21) with (S ee ) = (H hh ) — > —% 
for e — > ioo, corresponding to a constant density of states p(e) = 1/5 for |e| A. From 
(X), we find (Q) and hence the ensemble averaged supercurrent (/), 

oo N 

W = i kTA E E Si S n (^) ^ ( S " 6 (^)) - ^ (S e?l (^))] • (22) 

n=0 j=l 

Eqs. (21) and (22) contain all the information needed to determine the average supercurrent 
through a chaotic Josephson junction. 
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An analytic solution of Eq. (21) is possible in certain limiting cases. Here we discuss the 
case of high tunnel barriers, Tj <C 1 for all j. Then we may approximate Kj = (l/4)r., (A 2 — 
e 2 ) -1 / 2 and find 



(E e/l ) 
(Z he ) = Q*A 



(Z hh ) = -e(VA 2 - e 2 + E T ) 

n a 



\n\ 



n\ 2 A 2 -e 2 (VA 2 -e 2 + E T 



l A 2 -e 2 (VA 2 -e 2 + E T y 

1 



\Q\ 2 A 2 -e 2 (VAi-e2 + E T y 



N 



N 



J'=l 



2tt 



, (23a) 
(23b) 
(23c) 
(23d) 



j'=i 



The energy £"t is related to the mean dwell time through E^ = ft/2rd we ii [33]. The excitation 
gap in the spectrum of the Josephson junction is of order |fi(</>)| when r dwc n ^> K/A [18]. 
Substitution of Eq. (23) into Eq. (22) yields the supercurrent 



(/) = 2 -^kTGE T ±^L 



Q\ 2 A 2 + ut 



A 2 + ul + Eh 



where 



G = 



2e 2 E£i E 



■j=JVi+l 1 i 



(24) 



(25) 



is the conductance of the Josephson junction when the superconductors are in the normal 
state. 

For arbitrary transmission probabilities Tj, it is necessary to solve Eq. (21) numerically. 
We have studied the case that both point contacts have an equal number of modes (Aq = 
N 2 = N/2), and that all transmission probabilities are equal (1^- = T for all j). The average 
supercurrent at zero temperature for T = 0.1 and T = 1 is shown in Fig. 2. 



IV. ERGODIC REGIME 

The general result (21)-(22) describes the supercurrent in the ergodic regime r erg <C h/A. 
Within this regime, we can distinguish two further regimes, depending on whether the dwell 
time Tdweii = Ti/2Et is short or long compared to h/A. We discuss these two regimes in two 
separate subsections. 



A. Short dwell-time regime 

In the short dwell-time regime (when Td we ii h/A, or equivalently Et ^> A), the 
magnitude of the critical current I c = max^ I((f>) is set by the energy gap A in the bulk 
superconductor: I c ~ GA/e at zero temperature. The temperature dependence of I c can be 
neglected as long as kT <^ A, i.e. for temperatures T much less than the critical temperature 
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FIG. 2: Average supercurrent at zero temperature, computed from Eqs. (21) and (22) for the case 
N\ = N2 = N/2, Tj = r for all j. Left panels: T = 1; right panels: T = 0.1. The upper panels 
show (I) in the short dwell-time regime for E^/A = 1 (bottom curve), 10, and 100 (top curve). 
The bottom panels show (I) in the long dwell-time regime for E^/A = 0.01 (top curve), 0.1, and 
1 (bottom curve). The conductance G = (2e 2 /h)NT /4 and the Thouless energy = NTd/2-rr. 
Notice that (I) is in units of GA/e in the top panels, and in units of GE^/e in the bottom panels. 



T c of the bulk superconductor. In the case of tunneling contacts, evaluation of Eq. (24) with 
E T > A > kT yields 

= GA if( 7 sin(0/2))sin0 
e v /l- 7 2 s i I1 2 (0/2) ' 

The conductance G was denned in Eq. (25), the function K is the complete elliptic integral 
of the first kind, and we abbreviated 

(JVi n \ 1 / 2 /n \ _1 

E E r - r d E r * ■ < 27 > 
i=l j=N!+l / \k=l J 

The parameter 7 equals 1 for two identical point contacts with mode-independent tunnel 
probabilities. 

The result (26) could also have been obtained directly from the general formula for the 
zero-temperature supercurrent in the short dwell-time regime [19], 

eA f 1 , . . t sin0 , . 

{I)= ?h dtp{t) h , • \m m> (28) 

2ft Jo a/1 - tsm^(0/2) 

which relates (/) to an integral over the transmission eigenvalues t of the junction in the 
normal state, with density p(t). The transmission eigenvalue density for a chaotic cavity 
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with two identical tunneling contacts (Ni — N 2 — N/2, Tj = T j+N / 2 for j = 1, 2, ... , N/2) 
is given by [34] 

P (t) = y r i( 2 - r i) = , (29 ) 

One can check that the integral (28) equals Eq. (26) with 7 = 1 if T 3 1 for all j. For 
two identical ballistic point contacts (N\ — N 2 — N/2, Tj = 1 for all j), the density is 
p(t) = N(2n)- 1 [t(l - t)}' 1 / 2 [35, 36], which yields 

(I) = 2eA x G^fiarsinh[tan(0/2)Ucotan(0/2)Y (30) 
i7msin(0/2) V / 

Here G = N/A and E is the elliptic integral of the second kind. 
B. Long dwell-time regime 

In the long dwell-time regime (when 7d we n ^> h/A, or equivalently E? A), the magni- 
tude of the critical current is set by the Thouless energy, but retains a logarithmic dependence 
on A: I c ~ (GE T /e) ln(A/_E T ). The temperature dependence of J c can be neglected as long 
as kT E T . If kT 3> E-? (but still T <^ T c ) the critical current decreases, though only 
logarithmically: I c ~ (GE T /e) ln(A//cT). For the case of tunneling contacts, we find from 
Eq. (24) the expressions 

(/) = ^sin01n ( - 2A/Et ) kT « Br, (31a) 
e I Vl-7 2 sin 2 (0/2)i 



(j) = sm< 



111 (^) + CEulcr 



kT > E T , T < T c , (31b) 



where CEuier ~ 0.58 is Euler's constant. For ballistic contacts, we do not have a simple 
expression as Eq. (31), but the parametric dependence of / on A, E^, and kT is the same 
as for tunneling contacts (cf. Fig. 2). 

The logarithmic dependence on A of the supercurrent when E^ A arises because the 
Thouless energy E^ is not an effective cutoff for the Matsubara sum (8) or, equivalently, for 
the energy integration (7). Spectral correlations exist up to energies of order h/r CTg ^> E T . 
These long-range spectral correlations are responsible for the weak decay S eA oc \/u of the 
self-energy and p — S^ 1 oc 1/e 2 of the density of states. The superconducting energy gap A 
has to serve as a cutoff energy for the otherwise logarithmically divergent Eqs. (7) and (8), 
which explains the logarithm In A in Eq. (31). 



V. NON-ERGODIC REGIME 

When r crg > h/A, a random-matrix theory of the Josephson effect is no longer possible, 
because the non-ergodic dynamics on time scales shorter than r crg starts to play a role. 
To study the average supercurrent in this non-ergodic regime, we return to Eq. (8). On 
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substitution of Eqs. (1) and (2) we obtain an expression for I in terms of the scattering 
matrix S of the normal region, 



2vr °° 

/ = — kTj2^M, (32a) 

e n=0 

F(oo) = -^-4rte ln[l-a(«w) 2 %)e !$ 5*(-^)e- ,$ ]. (32b) 
h a<p 

The evaluation of the scattering matrix at the imaginary energy ioo n is equivalent to the 
evaluation of the scattering matrix at the Fermi level in the presence of absorption, with 
rate l/r abs = 2cu n /h = (2n + l)2irkT/h. We first consider temperatures kT ^> E T . Since 
uu n ^> Et = h/2Tdweii for all n in this high temperature regime, absorption is strong, r a b s <C 
r dwcii- The formal correspondence between Matsubara frequency and absorption rate helps 
to understand that, to lowest order in r a b s /Td W cii = E-p/u, the diagonal elements of S(icu) are 
given by the reflection amplitudes of the tunnel barriers in the contacts, Sjj = (1 — T^) 1 / 2 , 
while the off-diagonal elements satisfy 

■p "p poo 

(\Sij(iu)\ 2 ) = —jJ- / drP y (r)exp(-2a«-/R), i * j. (33) 

Z^fc=l i k JO 

The function P^ is the classical distribution of dwell times for particles that enter the 
quantum dot through mode j and exit through mode i. Because of the smallness of 
{\Sij(iuo)\ 2 ) = 0(Er T /uj), it is sufficient to keep only the lowest order term in an expan- 
sion of (F(cu)) in the off-diagonal scattering matrix elements, 



(F( \\ -—X^ V 2a{iu) 2 {\S{iu) ij \ 2 )sin0 

h " h ^^[l-aN^i-r^ii-aM^i-r^r 



Eqs. (32)-(34) permit a semiclassical calculation of the average supercurrent in the non- 
ergodic regime for temperatures kT ^> Ej<, where random-matrix theory fails. The only 
input required is the classical distribution of dwell times. 

On time scales greater than r crg , the distribution Py is exponential with the same mean 
dwell time Td we ii = h/2E T for all i, j: 

Pij(r) = exp(-2E T r/h). (35) 

The non-chaotic dynamics on time scales shorter than r crg enters through a non-universal 
form of Pij for r < r crg . We consider the case of a ballistic dynamics (size L of the normal 
region much less than the mean free path £). The ergodic time r crg ~ L/v-p is then a lower 
cutoff on P^, since the minimum dwell time L/vp is the time needed to cross the system 
ballistically. A qualitative estimate of (/) is obtained if we set Pij{r) = for r < L/v-p and 
approximate it by Eq. (35) for larger times. Substitution of this dwell-time distribution into 
Eq. (33) gives 

(\Sij(iu;)\ 2 ) = Jf J ^eM-^L/hvp), u » E T , 1 ± j. (36) 
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We next compute (F(cu)) from Eq. (34), replacing a(iuj) by its value —i for uj <C A. The 
result is 

G Et 

(F(uj)) = -exp(-2a;L/^ F )sin0, (37a) 

r _ 2e" g& 41^(2 - rp-H2 - r 3 -)-i 

Notice that G = G for the case of high tunnel barriers (Tj <C 1 for all j). We can now 
calculate the average supercurrent from Eq. (32). Eq. (37) is valid for Et C w « A, 
E T <C hvp/L <C A, and is sufficient to determine the supercurrent in the temperature range 
E T <^kT < A. Substitution of Eq. (37) into Eq. (32) gives 

G Et 

(I) = sine/. In (hvp/nkTL) , E T < kT < hv F I L < A, (38a) 

e 

(J) = — ^-2sin0 ex.p(-2irkTL/hv F ), E T < hv F / L < A;T < A. (38b) 

Eq. (38) has the same temperature dependence as the result of Ref. 11 for the double-barrier 
SNS junction. 

We now turn to low temperatures kT ^ E T . In this temperature regime, the Matsubara 
sum (32) contains terms with uj n ^ Et, for which the off-diagonal scattering matrix elements 
Sij{iu) n ) are not small and the approximation (34) is no longer valid. However, since Et <C 
h/r CTg , these Matsubara frequencies are well within the validity range of random-matrix 
theory. Therefore, we can use the results of Sec. Ill to compute (F(u)) for uj < Et and the 
semiclassical formula (34) for uj > Et- These two results match at uj ~ E T , because the 
validity range uj <C h/r CTg of random-matrix theory and the validity range uj Et of the 
semiclassical theory overlap (assuming r crg <C Td wc ii = K/2Er). 

For the case of high tunnel barriers, random matrix theory gives [cf. Eq. (24)] 

, „ . , . GEt sin 6 , , . , 

{FM)= M*)\>*S w<<R " F/i<<A ' (39) 

while the semiclassical formula (34) gives 

(F(uj)) = GEt Sm - ex.p(-2ujL/hv F ), E T <w<A. (40) 

UJ 

[The function £l(<f>) was defined in Eq. (23d).] The two results (39) and (40) have a common 
range of validity E T <C uj <C hv F /L, within which they can be matched. The result is a 
formula valid for all w< A, for a ballistic quantum dot with high tunnel barriers: 

^'^tSf^^- 2 *'' w<<a - (41) 

After substitution of Eq. (41) into Eq. (32) we obtain the average supercurrent in the low- 
temperature regime, 



GEt . 
(I) = sm< 



111 I / ^! 9 ==== ) - CEulcr 

LEt^I -7 2 sin 2 (0/2) 



kT < E T < hvp/L < A(42) 



[The parameter 7 was defined in Eq. (27).] The results (38) and (42) cover the entire 
temperature range below T c . 
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l c x e/G 


weak couplii 
ergodic (r erg <C ^/A) 


lg (Tdwcll > T crg ) 

non-ergodic (r crg » ^/A) 


strong coupling 

(Tdwell — T e rg) 


short dwell time 

(Tdwcll < h/A) 


A 




A 


long dwell time 

(Tdwcll > A) 









TABLE I: Parametric dependence of the zero-temperature critical current I c on the three time 
scales Tdwell, T erg , and h/A. The Thouless energy Ei = h/2r^ Ne \\. 

VI. CONCLUSION 

In Table I we summarize the parametric dependence of the critical current at zero tem- 
perature on the three time scales Tdweii, T erg , and h/A. We show the three new regimes for 
a weakly coupled normal region (T dwc n T crg ), and have included for comparison also the 
two old regimes for a strongly coupled normal region (T dw cii — T crg ). Apart from a logarith- 
mic factor, the critical current is given by I c ~ (G/e) mm(^,/rd we ii, A) in each of the five 
regimes. There is an additional logarithmic dependence on min(Td W eii/T erg , T& wc \\A/h) in two 
of the three new regimes. Upon raising the temperature, the critical current is suppressed 
at a characteristic temperature given by min(^/r erg , A). At lower temperatures, I c has a 
logarithmic T-dependence as long as T > h/r^ew and becomes T-independent at still lower 
T. 

In this work, we did not address the sample-to-sample fluctuations of the supercurrent, 
but calculated only the ensemble average. For strongly coupled diffusive Josephson junctions 
(rdwcii — T crg , L 3> I), the root-mean-squared of the fluctuations is a factor e 2 /hG smaller 
than the average critical current [19, 37]. Preliminary calculations in the ergodic regime 
indicate that the same is true for weakly coupled Josephson junctions (Td wc ii 3> T crg ), i.e. the 
r.m.s. fluctuations of I c are given by the entries in Table I, times e/h. 

We close with a remark on quantum dots with an integrable classical dynamics, such as 
rectangular or circular ballistic cavities. For energies e ^5 Et, the excitation spectrum of 
an integrable Josephson junction is quite different from its chaotic counterpart [18]: The 
density of states p{e) of a chaotic cavity in contact with a superconductor shows a gap of size 
E T around the Fermi level e — 0, while p{e) vanishes linearly when e — > for a rectangular 
or circular cavity. It is an interesting open problem, to compute the supercurrent through 
an integrable cavity and compare with the results for the chaotic case obtained in this paper. 
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